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Abstract 

We introduce the notion of hidden quantum correlations. We present the mean values of observables 
depending on one classical random variable described by the probability distribution in the form of 
correlation functions of two (three, etc.) random variables described by the corresponding joint pro¬ 
bability distributions. We develop analogous constructions for the density matrices of quantum states 
and quantum observables. We consider examples of four-dimensional Hilbert space corresponding to 
the “quantum roulette” and “quantum compass.” 

Keywords: entanglement, hidden quantum correlations, information and entropic inequalities, qudits, 
noncomposite systems. 

1 Introduction 

Quantum correlation phenomena, like the entanglement [T] present in composite systems, for example, 
in the system of several qubits, are known to play an important role in developing new quantum technolo¬ 
gies, including quantum computing [2j. Strong quantum correlations in two-qubit systems responsible 
for the violation of Bell inequalities w were checked experimentally W- In [7], it was suggested to 
extend the notion of entanglement in order to relate this phenomenon to correlation properties of single 
qudits. 

The new entropic inequalities reflecting the presence of correlations and analogous to the subadditivity 
and strong subadditivity conditions known for bipartite and tripartite systems [SHE] were found for 
noncomposite systems like single qudits or multilevel atoms pun]- Examples of qudits, including 
j = 3/2, were considered in this context in [[18]422| . and the results obtained show that the correlations in 
composite systems and the correlations in noncomposite systems can formally be considered as identical, 
using a common mathematical framework. 

The aim of this work is to develop the approach for describing both classical and quantum correlations 
in composite and noncomposite systems, using the same scheme based on the application of invertible 
maps of integer numbers s onto pairs (triples, etc.) of the integers (j. k ) employed in [16] . Employing these 
maps, we demonstrate that a single variable and its statistical properties, such as mean values, can be 
considered as the properties of several random variables described by the corresponding joint probability 

Cased on the talk by M.A. Man’ko at the 22nd Central European Workshop on Quantum Optics (6-10 July 2015, 
Warsaw, Poland). 


1 



distributions and given in terms of the correlation functions calculated for these several random variables. 
We show this property for both classical and quantum systems. 

This paper is organized as follows. 

In Sec. 2, we discuss the means and correlations in classical systems. In Sec. 3, we present examples of 
four- and eight-dimensional probability distributions. In Sec. 4, we study quantum states, and in Sec. 5 
we consider in detail the case of N = 4 along with entropic inequalities. We give our conclusions in 
Sec. 6. 

2 Means as Correlation Functions 

Our aim now is to consider correlations in a single qudit as correlations in artificial multiqudit systems. 
We start from classical states. 

Following mm we consider a set of nonnegative numbers Ai,A 2 ,... ,An which, in turn, provides 

A s 

a set of other nonnegative numbers 0 < p s = —- < 1 satisfying the normalization condition 

E,=i A j 

Ps = 1 - The numbers p s can be interpreted as the probability distributions of one random variable. 
Let us measure the observable F(s). For each value of the integer s = 1,2,...,iV, one obtains the 
result of the measurement F(s). Repeating the measurement L times, where L is a large enough integer, 
one obtains such statistical characteristic as the mean value of the measured observable 

N 

(F) = J2psF(s). ( 1 ) 

S=1 

We consider s as a random variable , F(s) as an observable , and p s as the probability distribution of one 
random variable, which we call the state. The other statistical characteristics described by the highest 
moments like, for example, variances 

N N 2 

(F 2 ) ~ ( F ) 2 = Y^PsF 2 (s) ~ (X>^00) , (2) 

S=1 8=1 


can also be obtained. 

Formally, one has two functions F(s) and p s defined on the set of integers s = 1, 2,..., N. 

We call the function p s the state and the function F(s) the observable due to the following reason. 

If one has the continuous variable x (a position of the particle with the Gaussian distribution P(x)), we 
call x the random variable and the distribution P{x) the state of the system. We extend this terminology 
to a discrete variable s and the distribution p s . 

The mean value ( F ) and variance (F 2 ) — (F) 2 are the functionals given by Eqs. (JT|) and (|2|) . In 
principle, one can formally define the functions F(s) and p s , as well as the functionals (Q]) and ([2]) . 
without the probabilistic interpretation of these objects. 

On the other hand, there exist functionals determined not by both functions F(s) and p s but by only 
one function p s . For example, Shannon entropy [23] associated with the probability distribution p s is 
given by the expression 

H = ~ Ps ln ps ; (3) 

S 
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the entropy being the functional of the state. 

Meanwhile, the entropy H is the functional which can also be considered formally without a pro¬ 
babilistic interpretation of the numbers p s . We point out the possibility to treat the functionals (F), 
(F 2 ) — ( F ) 2 , and H as objects that can be considered without their probabilistic interpretation because 
the numerical properties of these and other analogous functionals, e.g., all highest moments 

N 

(F k ) = Y^PsF k (s), (4) 

8=1 

like equalities and inequalities for these objects, exist independently of their relation to probabilities. 

One can repeat the above consideration for nonnegative numbers pj k ■ This simple observation can 
be used for obtaining some new equalities and inequalities for functionals (entropies, correlations, means, 
variances, and covariances) associated with tables of nonnegative numbers 0 < pj k < 1, j = 1,2, ... ,n, 
k = 1,2,... ,m, and N = nm, since the table can be considered as a joint probability distribution for 
two random variables. 

Within the framework of the interpretation of the table pjk as a joint probability distribution, the 
characteristics like entropy, mutual information, etc. naturally appear. These characteristics are known 
to satisfy the entropic inequalities for bipartite classical systems; see }2T |. 

On the other hand, the numerical expressions of these inequalities are valid independently of the 
probabilistic interpretation of the numbers in the table pjk- We employ this fact for obtaining new 
inequalities for the state (the probability distribution p s ) associated with one random variable and the 
function F(s) ( observable F(s)). The key tool to achieve this result is introducing the map of integers, 
namely, for s = 1,2,...,1V we construct the invertible map 

101,1; 2o2,l; ...; non, 1; n + lol,2; n + 2o2,2; ...; N — 1 o n — 1, m; JVon,ra. 

This map could be described as a procedure for introducing the function s(j, k). Such a function defined 
in the domain of integers j = 1,2,...,n and k = 1,2,... ,m provides for each pair of the integers j, k the 
value of the function equal to the integer s. The function is constructed using the invertibility condition; 
this means that for each value of the integer s one has only one pair of integers j, k corresponding to this 
value. Such construction was used in [16] to derive new entropic inequalities for qudit states. Here, we 
extend this construction to study the properties of observables associated with functions F(s). 

In fact, these observables associated with one random variable can be treated as observables connected 
with two random variables. To demonstrate this fact, we define the function <h(j, k) = F(s(j,k)). One 
can choose this function in the product form 

= </>{j)x(k). (5) 

The form of observable F(s(j,k )) provides the possibility to interpret the observable as the existence 
of two observables 4>(j) and x(fe) associated with two random variables j and k. Also the probability 
distribution p s can be chosen in the product form pj k = UjVk', this representation can be chosen with 
high ambiguity. In view of this representation, one can rewrite formula m for the mean value (F) as 
follows: 

N n m 

(F) = ^2 Ps F(s) = EE ${j> k )Pjk- (6) 

s=l j =1 k =1 
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Also one can introduce marginal probability distributions: 


n? = Y' p F = Y Ps UF): j = 1,2, (7) 

k= 1 k=l 

n n 

v k = YPjk = Y Ps U,k)’ k = 1 , 2 ,..., m. ( 8 ) 

j= 1 3=1 

If the numbers determining the joint probability distribution are such that pjk = IIj'P/., where 
Ilj = J2k=i = 1) one has for the marginal distributions © and (HI) the case of the absence 
of correlations between the observables associated with the function d>(j, k) given by ©. 

If the function F(s(j, k )) has the product form analogous to ([5|), the mean value of this function (F) 
given by ([6]) and written as 

n m 

( F ) = Y Y tti)x(k)p jk (9) 

3=1 k= 1 

can be interpreted as the correlation function, i.e., 

(F) = (cKj)x(k)). (10) 

Thus, for one random variable F(s) we obtain the formula for its mean value in the form of correla¬ 
tion function associated with two observables depending on random variables <p(j) and x(k), using the 
averaging procedure determined by the joint probability distribution pp-. 

In the case of integer N = ni^ns, where the factors in the product are integers, one can use the 
invertible map of the integers s onto the triples of integers j, k, 7, where j = 1, 2,..., m, k = 1,2,... , ri 2 , 
and £ = 1,2,..., n^. This means that we construct the function of three variables s(j. k, i) such that for 
each three integers we have only one integer s, and for each integer s we have only one triple of integers 
j, k, £. 

In view of this invertible map, the probability distribution p s used to describe statistical properties 
of observable F(s) depending on one random variable F(s) may be interpreted as the joint probability 
distribution p s (j^,i) = Pjke of three random variables. For this, we define the function T(j,k,£) = 
F(s(j,k,£)^, which can be chosen in the product form 

T(j,k,i ) = a(j)b(k)c(£). (11) 


Then one can write the equality 


N ni ri2 ri3 

(F) = Yps f {s) = LEE to\m)p s(j -,m) 

S=1 j = 1 fc=l t= 1 


( 12 ) 


or 


which means 


ni ri 2 n 3 

( F ) = YYY a (ti h ( k '> c ^)ps(j,k/)’ 

j =1 k =1 l=\ 

(F) = (a(j)b(k)c(£)). 


(13) 

(14) 
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Thus, we presented the mean value of the observable depending on one random variable in the form of a 
correlation function of observables depending on three random variables. 

Analogous representations can be developed for highest moments of the observable of one random 
variable. 

3 Examples of N = 4 and N = 8 

We recall that in our approach the integer s is the random variable , the numbers p s (the probability 
distributions) are the states, and the function F(s) is the observable, which has a value equal to the 
number F(s). For s = 1, 2,... , N, we have N values of random variable. One can use any other notation 
for the states and random variables, using an invertible map of the integers 1 , 2 ,..., iV onto another set 
of numbers mi, m 2 ,..., mjy. 


3.1 Case of N = 4 

We study the suggested construction on the example of N = 4. As an example, we consider these 
numbers as numbers associated with a casino roulette (or geographic compass). 

This means that we have four different positions of the casino roulette s = 1, 2, 3,4 (or four directions 
of the compass arrow). 

We use the map 1 -H- 1, 1; 2 -H- 1, 2; 3 «->• 2, 1; 4 -H- 2, 2 to label the four roulette positions by 
four pairs of numbers p jk (j, k = 1, 2), i.e., pi = pu, p 2 = p\ 2 , P3 = P 21 , and p 4 = p 22 . 

Now we introduce the observable F(s), which is a function of a random variable s equal to the number 
F(s) at each value of the variable. In this way, we have four numbers F(s = 1) = F(l), F(s = 2) = F( 2), 
F(s = 3) = F( 3), and F(s = 4) = F( 4). 

The mean value of the observable reads 

(F) = PiF(l) + p 2 F{ 2) + p s F(3) + p 4 F* (4). (15) 

The mean value (F) is a functional that depends on two functions p s and F{s), i.e., the state and 
observable. 

Using the mapping procedure developed, we can rewrite Eq. (fl5l) as follows: 

(F) = Pi\F(l, 1) + pi 2 F( 1, 2) + p 2 iF(2, 1) + p 22 F{2, 2) (16) 

or 

2 2 

= ( 17 ) 

j =1 k =1 

Now we choose the function F(j, k) in the form 

F( 1,1) = V»(l)x(l). ^(1, 2) = p(l)x(2), F( 2 , 1 ) = ^( 2 ) X ( 1 ), F( 2 , 2 ) = <^( 2 ) X ( 2 ). (18) 

One can introduce two other functions tp(j, k ) and x(j- k), which provide the same result of multiplication 

F(l,l) = ^(1,1) X (1,1), F(l,2) =^(l,2)x(l,2), F(2, 1) = <p(2, l)x(2,1), F(2, 2) = ip(2, 2) X (2,2). 

(19) 
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In fact, one should obtain these equalities if 


£(1,1) = ^(1), £(l,2) = p(l), £(2,1) = ^(2), £(2,2) = <p(2), 

X(1,1)=X(1), x(l,2)=x(l), X(2,1)=X(2), x(2,2)=x(2). 


( 20 ) 


We can interpret the functions <p(j,k) = £(s(j, fc)) and x(;j. k) = x{s(j,k)) as two specific observables 
or two different kinds of a function of one random variable. The results obtained can be summarized as 
the equality 

( F) = <£x>; (21) 

this means that (F), being the classical observable mean, can be interpreted as the correlation function 
of two classical observables £ and £• We call the correlations of these two observables £ and x depending 
on one random variable s the hidden correlations. 

Analogously, for N = n\n 2 • • - ng one can obtain the equality 


(F) = (<Pi<P2---<Pe)> ( 22 ) 

where the same (F) can be considered as the correlation function of i observables £i, £ 2 , - - -, £^ (hidden 
correlations). 


3.2 Case of N = 8 

Now we consider the case of N = 8, where we also have nonnegative numbers p±,p 2 , ■ ■ ■ Ps, with 
Ss=i Ps = 1- We use the map s -f-)- s(j. k, £), i.e., 

Pi='Pm, P2 = Fu2, P 3 = Vi2h Pa = F 122> Pb = Tbi 1, pe = V212, P 7 = F22i, Ps = ^222- 

The nonnegative numbers Vjkg satisfy the condition k t= 1 'Pjkl = 1- They can be interpreted as a joint 
probability distribution of three random variables j, k, and t. 

We turn to the observable F(s), s = 1, 2,..., 8. In terms of the probability distribution p s , the mean 
value (. F) reads (F) = F(s)p s . In view of the notation F[s(j,k,£)) = F(j,k,£), we arrive at 

(F) = J 2 V ^F(j,k,£). 

3 ,M 

If the observable F(j, k,£) is taken in the form 

F(j, k,£) = ip(j)x{k)u(£), 

we obtain 

( F ) = ^2<Pti)x(k)u(£)V jk g = (ip(j)x(k)u(£)). 
j,k,£ 

Thus, we obtain the result that the mean of a specific observable F(s) appears in the form of the 
correlation function of three observables 

<£(j,M) = p(j), x(j,k,£) = x(k), u{j,k ,£) = u{£), i.e., (F) = ( ££u). 
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The functions ip, Xi and u can be interpreted as observables depending on one random variable s. Thus, 
the mean value of the observable F(s) can be written as the correlation function of three observables 
£00, x00> and u(s), i.e., 

N N 

^ ~2psF{s ) = ^2p s p(s)x(s)u(s). 

S=1 s=l 

4 Quantum Qudit States and Observables 

In this section, we construct quantum states and observables, extending the approach discussed in 
the previous sections for classical systems. 

Given NxN matrix p ss i ( s,s' = 1,2,... ,7V). If p = p\ Tr p = 1, and p > 0, this matrix can be 
interpreted as the density matrix of qudit state with j = (N — l)/2. 

At N = nm, the matrix p ss i (s,s f = 1,2,... ,N) can be interpreted as the density matrix of two 
qudits with j\ = (n — l)/2 and j 2 = (m — l)/2, as well as at N = 711712713 , it can also be interpreted 
as the density matrix of three qudits with j\ = (n\ — l)/2, j 2 = {n 2 — l)/2, and j '3 = (713 — l)/2. An 
analogous interpretation can be provided for N = ni=i an d the matrix p ss t (s , s' = 1,2,..., N) can 
be considered as the density matrix of £ qudits with jk = (rik — l)/2. 

To provide such an interpretation, we use the map of matrix indices s ■<->■ j,k, s' ■<->■ i.e., 

s = s(j,k) and s' = s'(j',k'), while considering two qudits, and s = s(j,k,£) and s' = s'(j',k 'while 
considering three qudits, etc. We used this tool in [16]. In this paper, we study the possibility to extend 
this interpretation also for matrices of observables F ss / corresponding to the operators F acting in the 
Hilbert space %. 

We can write the matrices of observables either in the form 

Fss' Fs(j,k) s’ij' ,W) = Fjkj'k'i (23) 

or in the form 

Fss' = F s (j t k,£) s'(j',k',P) — Fjki^j'k'l'i (24) 

where indices j, k and j, k,£ take the same values as in the density matrix 

Pss' Ps(j,k)s'{j',k') = Pjk,j'k'i Pss' Ps(j,k,t) s'(j',k'— Pjkl,j'k'i'- 

Thus, both quantum states and quantum observables described by a density operator p and an observable 
operator F acting in the N x IV-dimensional Hilbert space H can be associated with the matrices p ss > and 
F ss 1 given in the basis | s), i.e., p ss > = (s \ p | s') or F ss i = (s | F \ s'). 

On the other hand, one can use the basis | s) =| s(j,k)) =\ j) \ k ), considering the Hilbert space H 
as the tensor product of two Hilbert spaces H = Hi (g) H 2 . In this basis, the matrix of the same density 
operator p reads 

Pjkj’k' = ( s(j,k) | p | s'(j',k'))-, 

this is the same numerical NxN matrix p ss r but with the matrix elements labeled by indices jk , j'k!. 
Analogously, for the observable F we can write the matrix (s \ F \ s') in the form 

(s{j,k) | F | s'(j',k')) = F jktj i k '. 
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Thus, we obtain the same NxN numerical matrix with matrix elements F ss > (s, s' = 1,2,..., N ) but the 
matrix elements are labeled by the indices jk and j'k! ( j,f = 1, 2 ,... ,n; k,k' = 1,2,..., m). The map 
introduced provides a chance to write the mean value of the observable F ss i = F s tj t k) s '(j',k') = Fjk,j'k' as 

N N n m n m 

(t ) = il Fp = ^2 ^2 Fss'Ps's = ^2 22 22 22 ^ s U-k) s'(j',k') Ps'(j',k') s(j,k) 

s=l s'=l j=l k=l j'=l k'=l 

n m n m 

— 22 ^2 ^22 Fjk.j’k' Pj'k',jk- (25) 

j =1 k =1 j '=1 k '=1 

If one takes the observable F in the form 


F = F 1 ^F 2 , (26) 

where F\ is the operator of the observable acting in the Hilbert space Tii and F 2 is the operator of the 
observable acting in the Hilbert space 'U 2 , Eq. (1251) reads 

n m n m 

h = EEEE ( F i)jj'{F 2 ) kk / Pj’k',jk■ (27) 

j =1 k =1 j '=1 k '=1 

In the case where F = F\ ® F 2 , we introduce two commuting observables 

-Ei = (A <g> l m ), F 2 (l n ®F 2 ). (28) 

For these two observables, the mean value of the observable F takes the form of the correlation function 
of the observables F\ and F 2 , i.e., 

(F) = (A-E 2 ). (29) 

As a result, we obtained a quantum analog of the classical probability relation (1211) . 

For N = nin 2 ■ ■ ■ ri£, we have 

(F) = {F\F 2 ■ ■ ■ Fg); (30) 

this relation is a generalization of Eq. (1291) . We showed that the same mean value (F) can be considered 
as the correlation function of £ commuting observables F p (p = 1, 2,... ,£). 


5 Entropic and Information Inequalities 


In this section, we consider the classical system with one random variable. We recall that there exist 
inequalities for entropies of joint probability distributions V(j, k ) of two random variables j and k of the 
form 


-22v(j,k) In V(j,k) 
jk 


< - 


22{[22 v ^' k ^ ]n [52 v ti> k ) ln l22 v (F k ) }. 


j k 


k j 






This inequality (the subadditivity condition) can be interpreted as the subadditivity condition for the 
probability distribution of one random variable 


- X p s lnp s < - X { [ X p s{j>fc) ] In [ X p am ] } ~ X { [ X 


, Ps{j,k) 


In 


3 


}■ 


The other example (for N = 4) is the possibility to consider the observable F(j. k ) for two random 
variables j and k, for example, F(l,l) = 1, F(l,2) = — 1, F( 2,1) = —1, and F( 2,2) = 1 as an 
observable for one random variable F(s) such as F( 1) = 1, F( 2) = —1, F( 3) = —1, and F( 4) = 1. In 
this case, 

(F) = piF(l) + p 2 F(2) + p 3 F( 3) + p 4 F(4) = p\ - p 2 - p 3 + Pa- 

Also 


(F) = V(l, 1)F(1, 1)+V{1, 2)F(1, 2)+V(2, 1)F(2, 1)+P(2,2)F(2, 2) = V{1,1)-P(1, 2)-P{2, l)+7>(2, 2). 
On the other hand, for two observables 

Fi(l,l) = l, Fi(l,2) = l, Fi(2,l) = -1, Fi(2,2) = -1 


and 

T 2 (l,l) = l, F 2 (l,2) = -1, F 2 ( 2,1) = 1, F 2 { 2,2) = -1, 
one has the correlation function of the form 

(FiF 2 ) = X V(j, k)A ( j , fe)F 2 (j, fc), 

j,k 

and this correlation function is equal to (F). In fact, 

4 

(FiF 2 ) = 1) - P(l, 2) - P(2, 1) + P(2,2)=p l -p 2 -p 3 +p A = Y J F(s)p s . (31) 

S=1 


In view of the invertibility of the applied map of the indices, we can introduce two observables T((s) 
and F.'(s), i.e., 

F((l) = l, F((2) = l, F((3) = -l, F[{ 4) = 1 

and 

F'(l) = l, F'(2) = -l, F'(3) = l, F'(4) = -l. 

Then one has the equality (F) = {F[F£) or 

4 4 

Y J F{s)p s = Y J F[{s)F! 2 {s)p s . (32) 


S=1 


1 = 1 


Thus, we showed that the mean value of observable F(s) can be interpreted as the correlation function 
of observables F\ (j. k) and F 2 (j. k). The subadditivity condition for functions p s of one random variable 
reflects the correlations of two artificial random variables (j, k) that are connected with two observables 
F\ and F 2 . Another interpretation of equality (132(1 reflects the fact that there exist hidden correlations 
of observables F[ (s) ad F^s) for the case of a single random variable s. 
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6 Conclusions 


To conclude, we point out the main results of this study. 

We showed that for a single qudit with j = (N — l)/2 it is possible to find commuting observables 
(Hermitian matrices), e.g., two observables p and x, such that the product of the observables provides 
the Hermitian matrix px = A. Then the mean value of the observable A can be interpreted as the 
correlation function of two observables {A) = (px) = Tr ( Ap ). 

In the case N = niri 2 ■ ■ -rik, one can find commuting observables p>i,p> 2 , ■ ■ ■ ,Pk, such that the means of 
the observable A = p\p 2 • ■ ■ Pk can be treated as the correlation function (A) = (pip >2 ■ • ■ Pk) = Tr (Ap). 

In the case of a multiqudit system with the same numerical density matrix p( 1,2,... ,k), the ob¬ 
servables pi, p> 2 , ■■■ ,pk have the physical meaning of the observables associated with each qudit in the 
composite system. 

Such an observation means that the quantum correlations known for observables associated with 
the subsystems are also available in single qudit systems like the quantum roulette and the quantum 
compass. We call these correlations the hidden correlations and they can be used in quantum technology 
applications analogously to the correlations in composite systems (like, e.g., the entanglement). 
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